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ABSTRACT 

This report presents an improved theoretical treat- 
ment of the two-frequency (24.5/73.6 MHz) rocket-to- 
ground CW propagation experiment for in situ electron 
density measurements. The propagation experiment measures 
the difference between two phase path integrals, and the 
local electron density is calculated from the rate of 
change of this difference. The basic relationships are 
derived directly from the concept of equivalent vertical 
electron content. This leads to formulas in which temporal 
and spatial ionospheric variations appear in explicit form. 
These variations give rise to errors, which depend upon 
both fluctuation rates and geometric factors associated 
with the rocket trajectory. Examples of these errors 
are given for typical experimental conditions. 


t 


INTRODUCTION 


Seddon’s two-frequency CW propagation experiment (Ref„ 1) 
has been used extensively for rocket measurements of iono- 
spheric electron density distributions. Basically, the 
experiment consists in phase path measurements accomplished 
by transmitting two CW signals from a moving rocket to 
receiving and recording stations on the ground. A suit- 
able low frequency was selected to obtain the ionospheric 
data and its 6th harmonic was used for comparison purposes. 

The experiment provides the difference between the high - 
and the low-frequency phase paths as a function of time (or 
equivalently altitude) from which the vertical electron 
density distribution can be computed. The 6-to-l frequency 
ratio was used for the experiments conducted from 1946 to 
1961. The lower of the two frequencies, initially 4,27 MHz, 
(1946) was changed to 7,75 MHz in 1954 and to 12,27 MHz in 
1959, Since 196? the experiment has been conducted primarily 
at 24,5 and 73.6 MHz using a 3-to-l frequency ratio. The 
24,5/73.6 MHz experiment has been performed on 8 rocket 
flights from Wallops Island, Va. (Ref. 2, 3 and 4), on 
6 rocket flights from a shipboard launcher near the coast 
of Peru (Ref, 5), and on 3 rocket flights from Thumba , India 
(to be published) , 

A general theory of the experiment which is applicable 
to the relatively low frequencies (4,27 and 7,75 MHz) used 
by Seddon has been described by Jackson (Ref, 6), At 24.5 
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and 73.6 MHz the formulas given by Jackson (Ref. 6) can 
be simplified as shown by Bauer and Jackson (Ref. 2). 

It was realized in previous treatments that the 
analysis yields an error term which arises from time 
variations in the ionosphere, and which affects the ac- 
curacy of the measurements primarily near the peak of the 
rocket trajectory. This has led to the conclusion that 
the measurements are usually reliable up to an altitude 
within 5-to-lO km of apogee. Below these altitudes the 
accuracy of the measurements was estimated to be about 5 
per cent. The present treatment re-examines this question 
in much greater detail than had previously been done. The 
present analysis starts'^ with basic principles but thkes 
immediate advantage of the simplifications which can be 
made at 24,5 MHz. This makes it possible to express the 
observations directly in terms of the electron density 
distribution, and to provide (at least to non-specialists) 
a clearer insight into the physical significance of the 
measurements. The present treatment retains the error 
term and provides an estimate of its magnitude as a function 
of rocket altitude. A 


Phase Path Concept 


A CW signal propagating from 
A to O (both points assumed to be 
stationary) will (if the path r = OA 





is in free space) experience a phase 
shift '4' given by: 


Fig. 1. Propagation path 
between two fixed points 
0 and A in free space. 


( 1 ) 


Y = — = r — cycles 
X ^ 

where : 

r = distance OA in km 
X = wave length in km 
f = frequency in cycles/sec. 
c = velocity of light in vacuo (km/sec) 
In the ionosphere the phase velocity Vp is: 



where : 

N 

n = refractive index = F(— g-, B, 0) 

f^ 

and , 
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N = electron density (electrons/cm ) , 

B = induction (gauss) of terrestrial magnetic field, 

0 = angle between B and direction of propagation. 

The parameter B causes the refractive index to be double- 
valued and this gives rise to two modes of propagation. The 
present discussion assumes that we are dealing with only one 
of these two modes. A more complete treatment will be given 
later in the Appendix. 

For the 24.5 MHz experiment, ray bending effects can be 
neglected and we can also assume that the direction of phase 
propagation coincides with the direction of energy propaga- 
tion. In other words we can assume that the wave is still 
propagating along AO, even if part of the path is within 
the ionosphere. The phase shift for an interval Ar is: 

f f 

ay = Ar — = — nAr cycles (2) 

P 


3 


( 3 ) 


and the total phase shift is therefore: 


f 

Y = — Jo n(r) dr cycles 


The integral n(r) dr, which is known as the phase path P, 
represents the distance r multiplied by ah average index n 
as seen by Eq. (4) : 


4. X 

p = Jo ^ “Jo 


CW Propagation Experiment (Basic Principles) 

The CW propagation experiment provides basically a 
measurement (as a function of time or equivalently r) of the 
difference between the high- and the low frequency phase 
paths, i.e: 


- '’l ■ I’d - J'o 

r 

■ J’o 


where the subscript H refers -to the higher frequency wave 
and the subscript L refers to the lower frequency wave. 

It should also be understood that n„ and n^ are both 
functions of r. For the frequencies undfer consideration 
(24.5 and 73 . 6 MHz) we can write (for N < 10® el/cc) 


n„ - n. = k N 
n Jj 


where k is a constant, 


k 


Although Eq. (6) is very nearly correct it is not used in 
the actual data analysis which takes into consideration 
both propagation modes and uses the exact relationship 

between phases indices and electron density. The actual ' 

procedure used to analyze the data is described in the | 

appendix where it is shown that the formulas derived in | 

the text require a minor modification when N is greater | 

than 3 x 10^ el/cc. j 

1 

Substituting Eq, (6) into Eq. (5) yields; | 

i 

I 

r r i 

P„ - P = f k N dr = k r N dr | 

H L o ‘Jo j 

j 

I 

The CW propagation experiment (at 24.5 MHz) can therefore 

r 

be viewed as yielding directly the quantity N dr, which ! 

is the total electron content from o to r (i.e, between 
points 0 and A of Fig, 1). The measured quantity U(r) 
will be defined for convenience as; 

U(r) - ^ CP„ - P,,) 

Thus 

r I'. 

N dr = U(r) (7) I 

where U(r) is the measurement provided by the CW propagation 

i 

experiment, ? 

1 

A knowledge of U(r) as a function of time makes it | 

i 

t 

possible to derive the electron density as a function of r j 
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(or equivalently the altitude z) as will be shown in the 
following sections. 

Theory of experiment for 9 - Constant 


Let point A in Fig. 1 represent the rocket -borne CW 
transmitter at a time t, and point 0 represents the re- 
ceiving site on the ground. We now let A move to a 
position as shown in Fig. 2a. We will assume that the 
motion is in the direction of OA, i.e. that the zenith 
angle 9 of the rocket position vector is constant. Let 
AA'^ = Ar and let At represent the time interval cor- 


responding to Ar, One example of 
the electron density N as a func- 
tion of r is given in Fig. 2b, 
When the rocket is at 
point A (OA = r) , the measure- 
ment indicated by Eq. (7) re- 
presents the cross-hatched area 
in Fig. 2b. Repeating the 
measurement when the rocket is 
at A'' (OA'^ = r + Ar) , and sub- 

tracting the previous measurement 
gives : 

U(r+Ar) - U(r) = area MQRS 

= Ar(Ng^) 

where N is the average density 

d V 

in the interval Ar, Thus 
" CU(r+Ar) - U(r) ]/Ar 

a V 



Fig. 2a. Geometry of experiment 
for 6 = constant. 



Fig. 2b. Electron density vs. r 


The above analysis (based upon a constant angle 0) also 
assumes that U(r) has not changed in value during the time 
interval At (corresponding to Ar) , i,e. it was assumed that 
U is a function of r only. The quantity as defined 
above represents also the average density for the height 
interval (z, z + Az) . The more general theory which does 
not assume that 9 is constant leads to a formula containing 
the expression: 

Cu(r+Ar) - U(r) ]/Ar , 

as well as other terms. It is therefore convenient to let 

n’ “ [U(r+Ar) - U(r)]/Ar (8) 

where N* represents the apparent density, l.e, the density 
corresponding to the assumptions: 

9 ■= constant; and U(r, t + At) = U(r, t) 

The assumption 9 = constant holds for most of the rocket 
ascent, because the rocket motion is then essentially radial 
with 'respect to the launch site (and with respect to the near-by 
receiving site) , The second assumption implies that the CW 
propagation experiment will yield most accurate results under 
quiet undisturbed ionospheric conditions. Ideal conditions 
for the measurements would be for the density to be a 
function of z only, i.e. independent of x (horizontal 
displacement) and t. Previous treatments of the subject 
have essentially assumed that the ideal conditions were 
fairly well approximated (during most of the ascent) by a 
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steep rocket trajectory, which permits the measurements 

to be made in a relatively short time (typically 5 minutes) 

and with a relatively small x displacement (typically less 

than 100 km) . The treatment for a horizontally stratified | 

ionosphere will be repeated here, using a more direct 

i 

mathematical procedure than was previously employed. I 

! 

I 

The concept of horizontal stratification is useful . ■ 

as a first order approximation to the actual situation for 
an undisturbed ionosphere. The general formula will there- 
fore include first order terms corresponding to horizontal 
stratification and second order terms arising from depar- 
tures from horizontal stratification. It will then be , 

shown that the second order terms have a negligible or 
very small effect during most of the rocket ascent, and 
that they can again have a fairly small effect over large 

i ' 

I 

portions of the descent trajectory, j 

I ’ 

In the subsequent discussion it will be convenient to i 

let : j 

r _ r ! ' 

U(r) = N(r) dr = N dr = N r (9) 

I ' 

where N is the average density from o to r, i,e, 

- 1 ' 

N = F J o 

We will now derive, for the case of horizontal stratification, > ■ 

i 

the following useful relationship; | " 

i 

X* z 

i Jo N(r) dr = 1 J^ N(z) dz (11) | 

i 

i 

I 
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where N(r) refers to the distribution along r and N(z) 
refers to the distribution along z. 

Horizontal stratification Implies that the electron density 
is only a function of the altitude z. Thus the density Nj 
can be considered to be a constant within the thin horizontal 
lamination shown In Fig. 3. 

Equation (10) can be written 


N = ^ 11m 


j=k 

E 


E 


Ar^o 

N » i lim 5^ N. Az ./cos 0 

^ j d 

N “ I Jo Ndz (12) 



Fig. 3. Constant density 
lamination in a horizontally 
stratified ionosphere. 


Theory of Experiment for Variable 6 and Horizontal Stratification 


The rocket is assumed to be 
travelling from Aj^ to Ag as in- 
dicated by Fig. 4, Let r^^ = OAj^ 
and rg “ OAg. The problem is to 
relate the density N in the height 
interval Zg - to the measure- 
ment U(rg) - U(rj^) . 



Fig. 4. Propagation paths 
for variable 0. 


9 


I 


Since we have assumed horizontal stratification we 


can v/rite for corresponding r and z: 


Z I* 

N(z) = i r Ndz = ^ [ Ndr 
z “ o r o 


Since U(r) = N r (see Eq. 9) we can write: 

= i- (Jo Ndz) r 

= _i_ Ndz 
COS0 o 

2 

U(rj) - U(rp - J ^ Ndz - ^ f 1 Ndz (13) 

2 o 1 o 

Using Eq. (8) to express AU and letting f = f f gives: 

o o 

^1 ^2 

= C— f Ndz + — r N dz 

COSGr, COS0T COS0O -i 

2 1 o 2 Zj^ 


or 


N'4i - ^ [■ 


_ 1 1 „ 1 „ Az 

At At '"cos0„ COS0, ^1 COS0O At" 


In the above equation, the expression inside the bracket 
represents the increment of the function l/cos0 as 0 increases 
from 0j^ to ^2 ^ ^1 ^ Letting At be the time interval 

corresponding to A0, it is seen that the above equation 
becomes for At 0: 


Nf ^1 dt ’008^^0=0^ + cos0j^ Zi 


( 14 ) 


_ r(+sin0)0T M , 1 XT i 

" z^ 

cos 0 0=0J^ 1 
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But z = r cos 6 

cos 8^ - rj^(sin6j^) 0^^ 

Nj^rj^ = rj^(tane^) + Nj^(r^ - r (tanQ^) 8^^) 

and more generally (dropping the subscripts 1) : 

N' = e N + N(1 - g) (15) 

where 

e = £t?2B8)i (16) 

r 

The approach used to derive Eq. (15) is based upon the 
geometrical interpretation of the measurement given by 
Fig. 4. A more formal procedure is to state that when 
At -^0, we can write: 


N- ^ ^ U(z,e,t) - ^ J-^ Ndz (17) 

Z 

Keeping in mind that J^Ndz does not depend upon 8 in 
this case, and expressing U in terras of partial deriva- 
tives gives: 


d ^ 5U dz du d8 dU dt 

dt 9z dt 90 dt 9t dt 

■ " t t " t t ' e t 

0,t=const z,t=const z,8=const 

The last term is zero for the horizontally stratified ionosphere 
(as defined earlier). Thus Eq. (17) becomes: 


« 
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N'r 


z — s- Ndz + e( Ndz) 

COS0 Sz O -Jo 


de 


+ z N 6 ^ 


cose 


's^ cose 


J. 

COS0 


which is the same as Eq. (14). 

A slightly more useful form of Eq. (15) is obtained 
as follows: 


N(l-.e) = N*(l-®) + ^ NV - ® N 
= N'(l-e) + e (N* - N) 
i.e. N = N’ + (N' - N) 


Equation (19) gives N in terms of N' , N and a factor 
e 

Tj — — which depends only upon trajectory parameters 

X 6 

(See Eq. 16) . 

It should be noted that N is for the altitude z 
(i.e. for point A^) and consequently it is known since the 
basic measurement provides N r^^. 

General Theory of the Experiment . 

The derivation of Eq. (19) assumed that the electron 
density was a function of altitude only. This assumption 
is not made in the more general treatment which will now 
be given. The basic measurement technique, as illustrated 
by Fig. 4, yields for the height interval Zg “ ^ 1 ’ 
quantity: 


2 

U(r 2 ) - U(rj^) = J N(r,t 2 , 02 )dr - J N(r,tj^, 0 ^) dr 


( 18 ) 


(19) 
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where : 

NCrjtgjQg) electron density distribution along OAg (for 

0 = ©2 and t = t^) , and 

N(r,tj^,0j^) is the electron density distribution along OAj^ 

(for 9 = and t = t^^) 

The density at a given altitude Z is not necessarily the 
same for the two propagation paths, because Gg / 0^^ and 
because t^ ¥ t^. The generalized expression for U must 
therefore include density variations due to t and 6, i.e. 

r 

U(r, t,0) = r N(r,t,0) dr 
' o 

The change of variable z = r cos0 will be used to 

r 

transform the expression N(r,t,9) dr, noting that t and 
0 are constants for the transformation. Thus 

dz = dr Cos e 

and 

r . z , , z 

N(r,t,9) dr = f N(z,t,0) = — L_ f N(z,t,e) dz (20) 

In Eq. (20) , the quantities t and 0 are functions of the 
upper limit of integration only. Furthermore N(z,t,9) 
does not represent the vertical electron density distribution 
along the Z axis. The function N(z,t,9) represents the 
distribution measured along a given path OA = r for the 
corresponding values of t and 0. 
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Thus : 


j* z 

U(r,t,0) = J^N(r,t,0) dr = J^N(z , t , 0)dz = U(z,t,0) 

or 

U(z,t,0) = N (z,t,0) (20a) 

Since N is now a function of 0 and t, the partial derivations 
of U in Eq. (17) with respect to 0 and t yield the following 
additional terms: 


_z ^ , z ^ 

COS0 d0 COS0 5t 


Letting ^ and recalling that -1-^ is multiplied by 0 

yields the following modified version of Eq. (15) : 


N’ = e N + N(l-e) + 




5Ns 

St'' 


SN SN 

where and are for constant z. 
0 0 ot 

Therefore Eq. (19) becomes: 


( 21 ) 


N = N' + 


1-e 


(N’-N) - 


(l-e)r 


(0 


Se 


at-' 


( 22 ) 


For very disturbed ionospheric conditions the error terms 
can be greater than N' , rendering the experimental measurement 
of N useless. For quiet and stable ionospheric conditions the 
error term is less than 5% for most of the ascent trajectory, 
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very large near apogee (where r is very small), and 5 to 

SN 

20% during descent. Some estimate of ^ for a quiet 
ionosphere can be derived from various measurements of 
the diurnal variation of N (Ref. 7 and 8). These 
measurements provide N for satellite-to-ground trans- 
mission paths and are therefore influenced primarily by 
the denser regions of the ionosphere at altitudes between 
200 and 500 km. Assuming that ^N/^t is proportional to 
N and independent of the altitude range, the above measure- 
ments (Ref. 7 and 8) indicate that 

J It ^ ^ 5U^ second) (23) 

Using the results of Little and Lawrence (Ref. 9), and 

making assumptions on 9N/S0 similar to those made on 

9N/St, one can arrive at the following upper estimate for 
dN 

I I < ^ N (per radian) (24) 

Thus l|^|< N and ||j| < Qg N 

For undisturbed ionospheric conditions and 

qg “ h‘ It is of interest to note that: 


(l-.e)r z 
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J * 

which follows from z = (rcosG) = r cos 0 - r(sin 0) 0 


dt 

z = (cos0) - r(tan6)0] 

z = r(cos0)^ - ^(tan6)0] 

r 

z = r(,cos9) (1-c) 

COS0 


r(l-e) z 


Thus the error term is Eq. (22) can be written: 


r If + If) I < I f -la - oi) 5 


The only difference in ascent and descent errors (at 
equal altitudes) is due to e, which is much greater for 
descent than for ascent. 

Limiting Cases . 

It is evident that Eq. (22) is not defined for either 
^ = 1 or 'r = 0. From the definition of e (Eq. 16) , it 
follows that e = 1 implies r = r(tan 0) 6 and from Eq. (26) 
it is seen that z = 0. Thus e = 1 corresponds to apogee 
when no measurement is made over a height Az* 

When r = 0 we cannot define N’ according to ^ = N*r. 

We can, however, write Eq. (14) as: 

r(tan0) § N - N r(tanG) 0 
= (N - N) r(tan0) 0 


( 26 ) 


(27) 
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which defines N in terms of N. The above expression may 
not be very helpful if r = 0 occurs very close to apogee, 
since z would then be very small, causing the error term 
to become dominant . 

SUMMARY 

The experiment yields as a function of time, or 
equivalently r(t) : 

r 

U(r) = N(r) dr = Nr (9) 

The electron density N at r(tj^) , i,e. for t = tj^ is given by: 

N = N' + (N’-N) (19) 

where 

, . £(ta^ ^ . 

r i 

N = - for t = t, 
r 1 

• • 

The quantities r, r, 9, 0,. tan 0, are derived from the known 
rocket trajectory. 

The error An in the measurement is: 

I^nI < I ^ (|e| qo + q.) N I (27) 

Z iS 0. 

where (for undisturbed ionospheric conditions) we can let: 

1 1 

qi = 5000 ^2 ^ ^ (9 in radians/sec.) 
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TYPICAL RESULTS 


The N(h) profile shown by Fig, 5 was obtained under 

undisturbed ionospheric conditions (K = 1) with NASA 8,25 

P 

(peak altitude: 629 km) using the 24.5/73.6 MHz experiment. 
The N(h) profile shown is based upon ascent data and it was 
calculated according to Eq. (19) , The large irregularities 
(AN/N) seen on the calculated topside profile are not con- 
sistent with typical topside sounder results which yield 
(particularly at midlatitudes) very smooth topside distri- 
butions, The AN/N fluctuations are consistent, however, 
with the errors indicated by Eq, (27) . Using in Eq. (27) 
the values of q^^ and qg given for undisturbed ionospheric 
conditions yield the following expression for the maximum 
per cent error AN/N: 


an _ z ^ N 


N 


5000^ N 


(28) 


The maximum error predicted by Eq. (28) is also shown (as a 
function of altitude) on Fig, 5. It is seen that the maxi- 
mum amplitudes of the fluctuations on the calculated profile 
are very consistent with the curve showing Thus one 

would conclude that the N(h) fluctuations are caused entirely 
by spatial and temporal variations along the propagation path. 

The N(h) profile shown by Fig. 6 was obtained during a 
disturbed period (Kp = 4) with NASA 18.98 (peak altitude: 

295.1 km) using the 24.5/73,6 MHz experiment. The profile 
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Fig. 5. Electron density profile and error profile for NASA 8.25 






(solid curve for ascent and dashed curve for descent) was 
calculated according to Eq. (19) . The calculations yielded 
a profile which was extremely smooth up to an altitude of 
260 km and very irregular above 260 km» The error analysis 
can again be used to determine whether or not the structure 
is real. Since the measurements were made under disturbed 
ionospheric conditions, the values of and q^ should be ' 

greater than the values used for the previous example. It 
was found that the computed error An/N would be comparable ; 

I 

I 

to the fluctuations seen on ascent above 290 km (i.e. near 
peak) if q^^ and q^ were both increased by a factor of 10. 

These increased values of q and q were therefore used to 

X a 

calculate the error curves (solid for ascent and dashed I 

i 

for descent) shown on Fig. 6. In spite of the large j 

increases in q^^ and the calculated error is signifi- 

i 

cantly smaller than the fluctuations shown at altitudes 
between 260 and 275 km. One would therefore conclude that 

{ 

the large AN/N values for altitudes between 260 and 275 km \ 

are real. It is also seen that below 260 km the shape 
of the descent profile is in very good agreement (except | 

that densities are slightly greater) with the shape of 
the ascent profile. The structure for altitudes between ’ ' 

t 

260 and 275 km is not the same on ascent and descent, and 

\ 

one would therefore conclude that this structure is localized. | 

This would contribute to the SN/d0 error and influence ; 

1 

measurements above 275. 

i 

J 

•i 

d 

■I 
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ALTITUDE (KM) 



10 ELECTRON DENSITY (EL/CC) 7^0 


Ascent and descent electron density profiles 
obtained with NASA 18.98, and corresponding 
error curves. The rocket was launched at 
0130 1ST (Indian Standard Time). 








APPENDIX 


Further Discussion of Measurement and Data Analysis Techniques 
The treatment of the experiment given earlier has neg- 
lected phase shifts introduced by the rotation of the trans- 
mitting antennas on the spinning vehicle. A relatively high 
roll rate of several rotations per second is normally used 
on sounding rockets to maintain their dynamic stability at 
altitudes above the drag region. This causes the plane of 
polarization of the transmitted wave to rotate at the rocket 
roll rate, and this rotation produces equal and opposite 
phase shifts on the two propagation modes. In order to 
cancel the phase shift due to rocket roll the basic phase 

measurement P„-Pt represented by Eq. (5) is conducted with 
n Jj 

both the ordinary and the extraordinary modes of propagation 
(Ref. 2). The sum of the two phase measurements (which is 
free of roll effects) is taken as the basic measurement. 

Thus the actual 24.5/73.6 MHz experiment (Ref. 3) yields 
as a function of time a net phase shift given by: 


^T ('•'jj ^^l^ORD. ^'“'l^EXT. 


where : 

H refers to the high frequency wave, 

L refers to the low frequency wave, 

ORD. refers to the ordinary mode of propagation, 

EXT. refers to the extraordinary mode of propagation, 
and the factor 3 is due to an effective frequency 
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multiplication by 3 of the received low frequency 
signal in order to compare its phase with the phase 
of the high frequency signal. 

Expressing in terms of the refractive indices gives: 

c~ >fo 
or 

’t ■ 5® J^O 

where Sn - (N„-NL)o„n + »*H-’'l>EXT 

For f = 73.6 MHz and f. =24.5 MHz, the quantity Sn is 
n ij 

essentially independent of 9 and B, and it is essentially 
a linear function of N, i.e. : 

Sn = k N 


Since k is not quite a constant, Sn will be written 


Sn = a k^ N 

where : 

k^ = 0.120, 

N is in units of 10® el/cc, and 
a is the function of aN shown in Fig. 7. 

For the typical N(h) profiles shown by Fig. 5 and 6, the 

0 

electron density is less than 10 el/cc, and the value of 
a would be: 1.0 < ^ < 1.04. Thus no serious error would be 

introduced by letting a = 1.0. 


(29) 


(30) 


24.5 MHz and 


Fig. 7. The function a for f^ = 
fjj = 73,6 MHz. ^ 


If desired a correction for a can be made as follows: 
Substituting Sn from Eq. (30) into Eq, (29) yields: 


7 = — k O' Ndr 

T c 0^0 


or 


r 

r a N dr 
o 


c7„ 




In the earlier treatment the quantity U was used to represent 
(c7^) /(f jjk^) . Retaining this definition of U, it is seen 


r 

U = PcvNdr = aNr 


that : 



and 


i- = ^ = (aN) apparent = a'N' (31) 

Thus the a correction consists simply in changing N to 
aN and it is readily seen that Eq. (19) would become: 


aN = q-’N’ + (a’N* - a N) 


where 

a’N’ 


1 ^ 

r " ^Vo 


(32) 


= (3 . 40x10^) Y^/r (Y in cycles/sec; r in km/sec) 

-- u 

^ “ "Vo 

= (3.40x10^) '•'rp/r ('J'rj, in cycles; r in km) 

The quantity a N given by Eq. (32) is converted to electron 
density using the appropriate value of a- The required value 
of a is found from Fig, 7 using the value of a N given by 
Eq. (32), 
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